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In this paper we report a measurement of the tt production cross section in pp colhsions at 
= 1.96 TeV using data corresponding to an integrated luminosity of 2.2 fb~^ collected with 
the CDF II detector at the Tevatron accelerator. We select events with significant missing trans- 
verse energy and high jet multiplicity. This measurement vetoes the presence of explicitly iden- 
tified electrons and muons, thus enhancing the tau contribution of tt decays. Signal events 
are discriminated from the background using a neural network, and heavy flavor jets are iden- 
tified by a secondary-vertex tagging algorithm. We measure a tt production cross section of 
7.99 ± 0.55 (stat) ± 0.76 (syst) ± 0.46 (lumi) pb, assuming a top mass nitop = 172.5 GeV/c^, 
in agreement with previous measurements and standard model predictions. 



I. INTRODUCTION 

In pp collisions at y/s = 1.96 TeV at the Tevatron, 
top quarks are produced mainly in pairs through quark- 
antiquark annihilation and gluon-gluon fusion processes. 
In the standard model (SM), the calculated cross section 
for top pair production at the Tevatron center-of-mass 
energy is 7.46t[5:^J^ pb [J for a top mass of 172.5 GeV/c^. 
This value can be enhanced by new processes beyond the 
SM such as top pair production via new massive reso- 
nances [2], while the comparison of the top pair produc- 
tion cross section in different decay channels can be sen- 



sitive to the presence of top decays via a charged Higgs 
boson 0. Thus, a precise measurement of the top pair 
production cross section is an important test of the SM. 
Both CDF and DO have performed many measurements 
of this quantity in different tt final states: the most recent 
published results, both measured in the decay channel 
with leptons and jets assuming mtop = 172.5 GeV/c^, 
are 7.70 ± 0.52 pb for CDF and 7.78 toll pb for 
DO g|. 

As the Cabibbo-Kobayashi-Maskawa matrix element 
Vtb is close to unity [6j [7] and the top mass mtop is 
larger than the sum of the W boson and bottom quark (6) 
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masses, in the SM the t — Wb decay is dominant and 
has a branching ratio of about 100%. Since the W sub- 
sequently decays either to a quark-antiquark pair or to 
a lepton-neutrino pair, the resulting top pair production 
final states can be classified by the number of energetic 
charged leptons and the number of jets. When only one 
W decays leptonically, the tt event is characterized by 
the presence of a charged lepton, missing energy due to 
the undetected neutrino, and four high energy jets, two 
of which originate from b quarks. In this lepton plus 
jets channel, one selects events with an energetic elec- 
tron or muon. For this paper we focus on an inclusive 
high-momentum neutrino signature of large missing en- 
ergy accompanied by jets. By not explicitly requiring 
leptons, our measurement is sensitive to all W leptonic 
decay modes including r decays of W^'s: about 40% of the 
signal sample obtained after the kinematic selection con- 
tains events with a t lepton in the final state. To ensure 
our measurement is statistically independent from other 
CDF results [5], we veto events with high-momentum 
electrons or muons as well as multijet events with no lep- 
tons {all-hadronic tt decays). This choice is expected to 
improve the final CDF combined cross section value: the 
previous 311 pb""'^ result [9j carried a weight of about 
17% in the combination [lO] . 

One of the major challenges of this measurement is the 
large background from QCD multijet processes and elec- 
troweak production of W bosons associated with b and 
c jets (heavy flavor jets), which dominates the signal by 
two orders of magnitude before any selection. In order to 
improve the signal to background ratio (S/B), a neural 
network is trained to identify the kinematic and topo- 
logical characteristics of SM tt events, and is applied to 
data to select a signal-rich sample. Top quarks are then 
identified by their distinctive decay into b quarks. Jets 
originating from b quarks (6 jets) are selected ("tagged") 
by their displaced vertex as defined by the SECVTX al- 
gorithm [TT]. After evaluating the average number of b- 
tagged jets for tt events using a Monte Carlo (MC) simu- 
lation, the number of signal events in the sample and the 
corresponding cross section are measured by counting the 
number of 6-tagged jets in the sample selected by the neu- 
ral network. The number of background ^-tagged jets is 
estimated using per-jet parametrized probabilities of &-jet 
identification, measured directly from data, rather than 
relying on theoretical prediction of cross sections and MC 
simulations. The results reported here are based on data 
taken between March 2002 and August 2007, correspond- 
ing to an integrated luminosity of 2.2 fb""'^, recorded by 
the CDF experiment at Fermilab. 

The organization of the paper is the following: Sec. |ll] 
contains a brief description of the CDF II detector and 
of the trigger requirements used for this analysis. The 
preliminary clean-up cuts applied to data are described 



are described in Sec. [V] and in Sec. VI respectively. We 
conclude the description of this measurement with a sum- 
mary of the different sources of systematic uncertainties 
in Sec. |VII[ while the cross section measurement is pre- 
sented in Sec. IVIIll 



II. THE CDF DETECTOR AND TRIGGER 
SYSTEM 

CDF II is a general-purpose, azimuthally and forward- 
backward symmetric detector located at the Tevatron 
pp collider at Fermilab. It consists of a charged-particle 
tracking system immersed in a 1.4 T magnetic field. The 
solenoid is surrounded by calorimeters and muon detec- 
tors [12]. The CDF II coordinate system uses 9 and 
as the polar and azimuthal angles respectively, defined 
with respect to the proton beam axis direction, z. The 
X-axis points towards the center of the accelerator while 
the y axis upwards from the beam. The pseudorapidity 
T] is defined as ry = — ln[tan(0/2)]. The transverse mo- 
mentum of a particle is px = psinO and its transverse 
energy Et = Es'mO. The missing transverse energy 
measures the transverse energy of the neutrinos via the 
imbalance of the energy detected in the calorimeters; it 
is defined by ~ \^t\ where = ^Jli^' 



in Sec. Ill followed by the discussion of the data-driven 
background parametrization in Sec. |IV| The kinematic 
variables characterizing the missing energy plus jets fi- 
nal state and the neural network based sample selection 



the 

index i runs over the calorimeter tower number and is 
a unit vector perpendicular to the beam axis and point- 
ing at the i-th calorimeter tower. The tracking system is 
composed of 8 layers of silicon microstrip detectors, ex- 
tending from 1.6 cm to 28 cm and covering up to I77I < 
2.0, surrounded by a 3.1 m long open cell drift chamber, 
providing |?7| coverage up to 1.0. Using information from 
the silicon detectors, the primary interaction vertex is 
reconstructed with a precision of ~ 15 /im in the plane 
transverse to the beam |13J . The energy of the particles 
traversing the detector is measured by electromagnetic 
and hadronic calorimeters segmented into projective tow- 
ers covering up to |ry| < 3.6. In the central region (|?7| < 
1.1) the calorimetric towers are 15° wide in and 0.1 in 77; 
in the forward region (1.1 < \ri\ < 3.6) the towers are 7.5° 
wide in azimuthal angle for I77I < 2.1 and 15° for \t]\ > 2.1. 
The electromagnetic section is made of lead-scintillator 
plates, while the hadronic section uses iron-scintillator 
ones. The transverse profile of electromagnetic showers 
is measured by proportional chambers and scintillating 
strip detectors. Muons are detected up to < 0.6 by 
drift chambers located outside the hadronic calorimeters, 
behind a 60 cm iron shield. Additional drift chambers 
and scintillator detectors provide muon detection up to 
I77I < 1.5. 

The CDF II trigger system [14] has a three level ar- 
chitecture designed to operate at 2.53 MHz and reduce 
the data rate to approximately 120 Hz to be written on 
tape. The data used in this measurement are collected 
with a purely calorimetric trigger, described in Sec. |II A| 
At Level- 1 (LI) calorimetric towers are merged in pairs 
along r] to define trigger towers: the LI can take a de- 
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cision on the energy contribution of individual trigger 
towers, on the sum of the energy of all the towers or on 
the missing transverse energy of the event. At Level-2 
(L2) trigger towers are merged into clusters by a sim- 
ple clustering algorithm ^15j, while at Level-3 (L3) jets 
of particles are reconstructed by a fixed cone-based al- 
gorithm |16j . the radius of the cone in the rj — (p plane 
being 0.4, 0.7 or 1 depending on 

the specific trigger. 

A. The Multijet trigger 

The data used in this analysis are collected by a multi- 
jet trigger. This trigger requires at LI the presence of at 
least one central trigger tower with Et > 10 GeV, and at 
L2 at least four calorimetric clusters with Et > 15 GeV 
each and a total transverse energy greater than 175 GeV. 
The latter threshold was 125 GeV before February 2005 
and was increased to reduce the trigger rate at higher 
instantaneous luminosity. Finally, at L3, at least four 
jets with Et >10 GeV {AR = 0.4) are required. This 
trigger was specifically designed to collect all-hadronic 
tt events, where the final state nominally consists of six 
jets, but has a large acceptance also on final states char- 
acterized by ^T ^iid high jet multiplicity. Moreover the 
collected sample is uncorrelated with those used for top 
cross section measurements in the lepton plus jets final 
state, which are selected by requiring the presence of a 
high momentum lepton. The choice of this trigger is also 
driven by the analysis strategy: the top cross section 
measurement is performed by counting the number of b- 
tagged jets from top decay in the final sample and the 
multijet trigger provides a sample that is unbiased with 
respect to the 6-tagging algorithm, as it does not apply 
any requirement on tracks. 

III. PRELIMINARY REQUIREMENTS 

Events satisfying the trigger requirements are used in 
the analysis only if they were collected with fully oper- 
ational tracking detectors, calorimeters and muon sys- 
tems, and if their primary vertex is located within ± 60 
cm along z from the center of CDF II detector. Jets 
and ^T corrected [17j for multiple pp interactions in 
the event, non uniformities in the calorimeter response 
along r], and any non-linearity and energy loss in the 
uninstrumented regions of the calorimeters. The ^t is 
also corrected for the presence of high-p^ muons. We 
consider jets with AR = 0.4, corrected transverse energy 
Et > 15 GeV and \ri\ < 2.0; the total transverse energy 
of the event ^ Et is defined as the sum of all jets Et- 

Events are required to have at least three jets and no 
central high-pT {pt > 20 GeV/c) reconstructed electrons 
or muons to avoid overlaps with other top cross section 
measurements in lepton and jets final states [18] . In the 
same way, overlaps with top all-hadronic analyses are 



avoided by rejecting events with low ^t significance ^t^, 
defined as ^ = ^t I where and ^ Et are 
measured in GeV, as the resolution on the ^t observed 
to degrade as a function of the total transverse energy 
of the event. The ^^'^ is typically low when ]fiT arises 
from mismeasurements, so events are required to have 
y^T^'^ > 3 GeVi/2. 

Throughout the paper, the impact of analysis require- 
ments on signal is evaluated using inclusive tt samples 
generated with pythia version v6.216 TH] and processed 
through CDF II detector and full trigger simulation 20] , 
assuming a top mass nitop = 172.5 GeV/c^ (the most 
recent CDF and DO combined result on top mass is 
rritop = 173.3 ± 1.1 GeV/c^ [H]). After the preliminary 
cuts, there are 94 217 events remaining in the sample 
with at least four jets, with an expected signal to back- 
ground ratio S/B of 1.4% and 44 310 events in the sample 
with exactly three jets, with an expected S/B lower than 
0.1%. The latter sample, with a very low tt signal con- 
tent, will be used to derive the background parametriza- 
tion in Sec. HVl 



IV. BACKGROUND PARAMETRIZATION 

In this analysis, top quarks are identified by their de- 
cays into b quarks. However, many processes other than 
the decay of the top quark can give rise to b jets, there- 
fore a procedure to estimate the number of 6-tagged jets 
yielded by background processes is needed. 

The background-prediction method used in this analy- 
sis rests on the assumption that the probabilities of tag- 
ging a b jet in tt signal and background processes are dif- 
ferent: the differences are due to the distinctive proper- 
ties of b jets produced by the top quark decays, compared 
to 6 jets arising from QCD and W boson plus heavy flavor 
production processes. In this hypothesis, parametrizing 
the &-tag rates as a function of jet variables in events 
depleted of signal allows one to predict the number of 6- 
tagged jets from background processes in any given sam- 
ple. The parametrization is derived from a sample of 
pure background: the data sample after the prerequisites 
requirement with exactly three jets, where the ti signal 
fraction is lower than 0.1%. In this sample the per-jet b- 
tagging probability is found to depend mainly on jet Et, 
the number of good quality tracks contained in the jet 
cone Ntrk, and the ^t projection along the jet direction 

^T^''^, defined as ^t^''^ = cos A(/)(^j,, jei). These 
variables are chosen for the parametrization and their 
corresponding tag rate dependence is shown in Fig.[l] We 
expect the tagging probability to depend on jet Et and 
Ntrk due to the implementation details of the 6-tagging 
algorithm |TT]. In more detail, the 5-tagging probabil- 
ity decreases at high jet Et due to the declining yield of 
tracks passing the quality cuts required by the 5-tagging 
algorithm; it increases if a greater number of good quality 
tracks is associated to the jet cone. The ^t projection 
along the jet direction is correlated with the heavy flavor 
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component of the sample and the geometric properties of 
the event. Neutrinos from the 6-quark semi-leptonic de- 
cays force to be ahgned with the jet direction, while 
neutrinos from W boson decays are more likely to be 
away from jets. For this reason the 6-tag rate is enhanced 
at high positive values of l^jP'^-' . 

A 3-dimensional 5-tagging matrix V is defined using 
the per-jet 6-tagging probabilities, and its binning is de- 
fined to avoid any entries with empty denominators. The 
matrix assigns the probability that a jet is 6-tagged given 
its Et, Ntrk and ^j.^"'' . 

The total number of expected background 6-tagged jets 
Nexp in a given data sample can be calculated by sum- 
ming the 6-tagging probabilities over all jets in the se- 
lected events: 



N 



N — 



E 

1=1 



jets,, 

E 



(1) 



where the index k runs over the number of jets Nj^ts in 
the i-th event and the index i runs over all the N^vents 
events in the sample. Due to the method we used, N^.^^ 
is obtained under the same assumption that allowed the 
parametrization of the tagging rate, i.e. a negligible tt 
signal content. For jet multiplicities greater than three, 
Nexp will overestimate the number of 6-tagged jets from 
background events, Nhkg, due to the presence of tt events 
in the sample. If Nobs is the number of 6-tagged jets ob- 
served in the data sample and we assume that the differ- 
ence between TVob^ and Ni,kg is due to the presence of tt 
signal, the number of top events in the sample ntop can be 
evaluated as ntop = {Nobs - Nbkg)/e^ag, where e^^| is the 
average number of 5-tagged jets per top event, defined 
as the ratio between the number of 6-tagged jets and the 
number of ti events and calculated from MC simulation. 
If Ugyt is the number of events in the sample, we can esti- 
mate the number of 6-tagged jets due to the background 
Nbkg, by rescaling the value N^xp predicted by the 6- 
tagging matrix as follows: Nbkg = Nexpinevt-ntop)/nevt- 
Putting the expressions for Nbkg and ntop together, we 
can evaluate Nbkg as the limit of the following iterative 
formula: 



Nbkg = lim 



(2) 



with 



N 



vt ~ {Nobs — Nf,xp,i-l)/£ 



exp,0 



tag 



where N, 



exp,Q 



Nf-xp is the number 



(3) 



fixed dur- 



ing the iteration - of expected 6-tagged jets coming 
from the tag rate parametrization before any correction. 
In our calculations the iterative procedure stops when 



< 1%. This correction procedure is used 
in all samples with more than three jets to remove the 



tt contribution from the background estimation. This 
method does not require knowledge of the top produc- 
tion cross section, once we assume that the dependence 
of £tag from the top quark mass is negligible. 

To ensure the correct behavior of the 6-tagging 
parametrization, an important check consists in calculat- 
ing the predicted number of 6-tagged jets in data samples 
with jet multiplicities higher than three, and compare it 
with the actual observed number of 6-tagged jets. The re- 
sult of this check is shown in Fig. [2j Taking into account 
the expected contribution to the number of observed 6- 
tagged jets due to the presence of tt events in the sam- 
ple {mtop = 172.5 GeV/c^), the agreement between the 
amount of observed and predicted 6-tagged jets is good 
in all the jet multiplicity bins, being exactly the same by 
definition for 3-jet events, on which the parametrization 
is calculated. 



V. SIGNAL AND BACKGROUND 
CHARACTERIZATION 

In the data sample under investigation, tt pairs are 
overwhelmed by multijet QCD and W + jets events. The 
main feature of the tt decay channel analyzed here is a 
considerable amount of ^^n, the only observable signa- 
ture of the presence of neutrinos from W leptonic de- 
cays. However, missing transverse energy can be also pro- 
duced by jet energy mismeasurements, by semi-leptonic 
decays of 6 quarks in QCD events and by the leptonic 
decay of the ly in -I- jets events. To discriminate 
against the possible sources of missing transverse en- 
ergy on a geometric basis we use the quantity A0„ii„ = 
min A(j){^rp, jet), defined as the minimum angular differ- 
ence between the and each jet in the event. This 
quantity is expected to be large in W ^ Iv decays and 
in tt ]^rp + jets events. On the other hand, for QCD 
background events where the main source of is due 
to jet energy mismeasurement, the l^rp is expected to be 
aligned with the jet direction and the value of A(/)„im 
close to zero. 

Other kinematic variables related to the topology of 
the event characterize the ti production with respect to 
background processes. Let Qj(j=i,3) be the eigenvalues of 

where a, 6 run over the three space coordinates, and Pj 
is the momentum of the jet j. The sphericity S is de- 



the normalized momentum tensor M, 



ab 



fined as S = ^{1 — Q3). S is zero in the limiting case 
of a pair of back-to-back jets, while approaches one for 
events with a perfectly isotropic jet momenta distribu- 
tion. The aplanarity A is defined as A = |(5i and lies in 
the range [0, Extremal values of A are reached in the 
case of two opposite jets and in the case of evenly dis- 
tributed jets, respectively. Jets emerging from a ti pair 
are expected to be uniformly distributed and, as a con- 
sequence, they will hardly lie on the same plane: thus we 
expect high aplanarity and sphericity values for events. 
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FIG. 1. b tagging rates as a function of jet Et (a) Ntrk (b) and ^^j,*""^ (c) for the data sample with exactly three jets in the 
event. 
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FIG. 2. Tagging matrix check on data after preliminary 
requirements and before any additional kinematic selection. 
Observed and predicted number of fe-tagged jets as a function 
of the jet multiplicity are shown in the figure, statistical er- 
rors only. The expected contribution coming from ti events 
based on the theoretical cross section of au = 7.45 pb is also 
shown. Background error bands are centered on the inclusive 
top plus background prediction. 



In addition to kinematic variables describing the topol- 
ogy of the event, distributions of energy- related variables 
can be useful to discriminate ti events over their back- 

V Et 

ground. The centrality C is defined as C = , vifhere 

\/s is the invariant mass of the jets system. In the case 
of tt pairs decaying hadronically, jets are emitted prefer- 
ably in the transverse plane {r — (j) plane) , so we expect 
to have a greater amount of energy emitted in this plane 
giving values of C closer to 1 with respect to background 
events. The variable ^3 Et is defined as the sum of all 
jets Et in the event except the two leading ones. In 
QCD events, the two most energetic jets are produced 
by qq processes, while the least energetic ones come from 
gluon bremsstrahlung; on the contrary, in tt events up 
to six jets can be produced by hard processes, and as a 



consequence ^3 Et can help discriminating signal and 
background. All these variables will be used in the fol- 
lowing section to train a neural network to discriminate 
tt signal from background. 



VI. NEURAL NETWORK BASED EVENT 
SELECTION 

In order to enhance the signal to background ratio in 
the data sample, we use a neural network (NN), trained 
to discriminate tt — )■ ^t + j^ts signal events from back- 
ground. The NN is built using the neural network imple- 
mentation in ROOT [23] . 

We apply an additional kinematic requirement on data 
with at least four jets, removing events with low angle be- 
tween jets and ^t with the cut A(j)„ii„ (^TT 3^^^) > 
In the data these events come mostly from mismeasured 
jets and are difficult to simulate in MC. They are re- 
moved from the NN training to prevent it from converg- 
ing on artificial differences between MC and data rather 
than real physics effects. After this requirement, we are 
left with 20 043 events in the sample with at least four 
jets, with an expected S/B of 3.5%: since this sample 
has low tt signal fraction and no overlap with data used 
to determine the background parametrization, it will be 
used as background training sample. As signal training 
sample we use the same amount of events passing the 
same event selection of data, taken randomly from the tt 
simulation. We use as inputs for the network the follow- 
ing kinematic variables, normalized with respect to their 
maximum values: the transverse energy of the leading 
jet E^, A(j)^,ni^TJet), ^t"' , E^t, Es^t, the cen- 
trality C, and the topology-related variables sphericity 
S and aplanarity A. The distributions of the input vari- 
ables used in the NN for both the signal and background 
training samples are shown in Fig. [3j 

After the training, the ^-tagged data are processed 
by the NN: Fig. |4] shows the number of observed b- 
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FIG. 3. Distribution of neural netw ork i nput variables for m ultijet da ta 
transverse energy of the leading jet E]^ 



topology-related variables sphericity S (g) 



A(t> 

and aplanarity A (h) 



(ba ckgro und) a nd t t MC (signal) s amp les: the 



jet) (b) ^j,"^ (c) J^^T (d) J^s^T (e) the centrality C (f) 



and the 



tagged jets versus the NN output NNo„t, along with the 
corresponding background prediction from the tag rate 
parametrization and expected contribution from tt sig- 
nal {rritop — 172.5 GeV/c^), for events with at least four 
jets and with exactly three, four and five jets. The good 
agreement between data and the sum of expected back- 
ground and tt induced 6-tagged jets in the high NN out- 
put region is both a confirmation of the effectiveness of 
the method we use to estimate the background and a 
demonstration of proper NN training and performance. 

In order to select a signal-rich sample to perform the 
cross section measurement, we choose to cut on the NN 
output value NNo^t- The cut is chosen with the aim of 
minimizing the statistical uncertainty on the cross sec- 
tion measurement, maximizing S/B where S is the ex- 
pected number of 6-tagged jets from the signal and B is 
the predicted number of background 6-tagged jets. The 



former quantity is evaluated from an inclusive ti MC sam- 
ple, while the latter is derived using the 6-tagging matrix 
parametrization on data. The result of this optimization 
procedure is a cut on NNo„t > 0.8, which gives the low- 
est expected statistical uncertainty on the cross section 
(~ 8%) and an expected S/B ^ 4. The expected signal 
sample composition after this cut is shown in Table |l] 



VII. SYSTEMATIC UNCERTAINTIES 

The top quark pair production cross section is mea- 
sured as: 

a{pp ^ tt) X BR{ti ^^^+ jets) ^ ^"'^ ~ .^^^^ (4) 
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FIG. 4. Observed and matrix-predicted number of b-tagged jets versus neural network output in the multijet data for all events 
with at least four jets (a) and for events with exactly three (b) four (c) and five [(d)] jets, along with the expected contribution 
due to tt signal events. Background error bands are centered on the inclusive top plus background prediction. 



TABLE I. Expected signal sample composition after the NNout > 0.8 cut. 



Number of jets 3 4 5 6 7 8 9 Total 

all-hadronic (%) 0.1 0.5 1.7 4.9 7.8 10.2 9.9 2:3" 

e +jets (%) 26.7 25.2 35.5 36.2 35.1 33.6 33.8 32.1 

/i +jets (%) 32.2 32.5 19.1 15.8 14.5 16.1 12.9 22.7 

dileptonic (%) 6.5 2.3 1.1 0.7 0.5 0.6 0.00 1.4 

hadronic r +jets (%) 15.8 21.9 30.5 31.5 31.8 29.9 34.4 27.7 

leptonic T +jets (%) 11.7 13.2 9.9 9.1 8.9 8.1 8.4 10.8 

TT (%) 1.1 1.1 0.6 0.5 0.4 0.4 0.3 0.8 

e/fi + T (%) 6.0 3.4 1.6 1.3 0.9 1.1 0.3 2.2 



where Nobs and N^xp are the number of observed and 
predicted tagged jets from background in the selected 
sample, respectively; ekin is the kinematic efficiency of 
trigger, preliminary requirements and neural network se- 
lection determined using inclusive MC tt events; ef^g, de- 
fined as the ratio of the number of ^-tagged jets to the 



number of events in the inclusive ti Monte Carlo sample, 
gives the average number of 6-tagged jets per tt event. Fi- 
nally, L is the integrated luminosity of the dataset used. 
All quantities in the denominator of Eq. |4j as well as the 
number of expected ^-tagged jets, are subject to different 
sources of systematic uncertainties. 
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The kinematic efficiency ekin is evaluated on inclusive 
ti samples generated with pythia and its uncertainty 
arises from the particular choice of the MC generator, 
the set of parton density functions used in the generator, 
as well as the modeling of color reconnection effects and 
of the initial and final state radiation. 

The MC generators differ in their hadronization 
schemes and in their description of the underlying event 
and multiple interactions: in order to evaluate the gener- 
ator dependence of the kinematic efficiency, we compare 
the value of ekin obtained using pythia with the value 
obtained on a sample of events generated with herwig 
v6.510 and take the relative difference as the sys- 
tematic uncertainty. 

The choice of parton distribution function (PDF) affects 
the kinematics of tt events, and thus the acceptance for 
signal events. We estimate this uncertainty comparing 
the ekin value derived from MC samples based on the 
default PDF CTEQ5L |25j with the one obtained using 
samples based on MRST72 and MRST75 dH, which dif- 
fer by the value of the strong coupling constant as used to 
compute the PDF; we also consider the difference in the 
value of ekin obtained with the leading order (LO) and 
next to leading order (NLO) calculations of PDFs, evalu- 
ated using defauh CTEQ5L (LO) and CTEQ6M (NLO) 
PDFs, and derive the corresponding uncertainty. We add 
these two contributions in quadrature to obtain the sys- 
tematic uncertainty due to our choice of PDF. 
Uncertainties arising from the modeling of color recon- 
nections effects are estimated by evaluating the shift in 
the kinematic efficiency using two samples of events gen- 
erated by PYTHIA, corresponding to different models of 
color reconnection 27 . 

Additional jets coming from initial and final state ra- 
diation (ISR and FSR), might change the sample com- 
position and affect the efficiency of the kinematic selec- 
tion. The systematic uncertainty related to these effects 
is evaluated by calculating the shift in ekin using inclu- 
sive ti samples with different amounts of initial and final 
state radiation. 

The systematic uncertainty due to the calorimeter re- 
sponse is accounted for by varying the corrected jet en- 
ergies within ±1(7 of their corresponding systematic un- 
certainty, and recalculating ekin after these variations. 
Finally, ekin is also affected by the simulation of the trig- 
ger requirements on MC events and a trigger acceptance 
uncertainty is determined by comparing trigger turn-on 
curves between MC and data events. 

The average number of 6-tagged jets per tt event is 
affected by the uncertainty on the scale factor used to 
account for the different efficiency of the 6-tagging algo- 
rithm in data and in MC. The systematic uncertainty on 
ef^g is obtained varying the scale factor within the ±la 
range from its central value and determining, on the MC 
sample, the difference in terms of average number of b- 
tagged jets per event with respect to the standard ef^g 
value. 

The uncertainty on the number N^xp of matrix-predicted 



TABLE II. Summary of relative systematic uncertainties on 
the signal efficiency, and other uncertainties related to the 
cross section evaluation. 



Source Uncertainty 


ffcin systematics 




Generator dependence 


3.9 % 


PDFs 


1.2 % 


ISR/FSR 


2.7 % 


Color Reconnection 


4.3 % 


Jet Energy Scale 


4.2 % 


Trigger simulation 


3.0 % 


Other systematics 




^tag (^tag scale factor) 


3.9 % 


Background parametrization 


2.5 % 


Luminosity measurement 


5.8 % 



TABLE III. Input values for the cross section measurement. 
e-tag is the average number of 6-tagged jets per top event. 



Variable 



Symbol Value 



Integrated Luminosity (p5"^) L 2207.5 ± 128 
Observed 6-tagged jets Nobs 636 
Background 6-tagged jets Ni,kg 131 ± 9.6 
Kinematic efficiency (%) ekin 3.53 ± 0.29 
6-tagged jets per tt event (%) e^^^ 81.1 ± 3.2 



^-tagged jets is calculated by comparing the number of 
^-tagged jets yielded by the tagging matrix to the actual 
number of observed SECVTX tagged jets in a control sam- 
ple depleted of signal (events with NN output lower than 
0.6). The relative difference between the expected and 
observed number of tagged jets is taken as an estimate 
of the uncertainty of our background prediction. 
The luminosity measurement is affected by two sources 
of uncertainty: the acceptance of the luminosity monitor 
and the total inelastic pp cross section (60.7 ± 2.4 mb). 
The uncertainties on these quantities are 4.2% and 
4.0% respectively, giving a total uncertainty of 5.8% on 
the integrated luminosity calculated for any given CDF 
dataset. 

The summary of all sources of systematic uncertainties 
to the cross section evaluation is listed in Table HIl 



VIII. CROSS SECTION MEASUREMENT 

After the neural network selection we are left with 1420 
events with at least four jets of which 636 are 6-tagged. 
Inserting in Eq. [4] the input parameters quoted in Ta- 
ble IIIIl the measured cross section value is 



ati = 7.99 ± 0.55 (stat) ± 0.76 (syst) ± 0.46 (lumi) pb 
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FIG. 5. Observed and matrix-predicted number of 6-tagged 
jets by jet multiplicity in the multijet data after cut on neu- 
ral network output NNout greater than 0.8, along with the 
expected contribution of inclusive tt signal normalized to the 
measured cross section. Background error bands are centered 
on the inclusive top plus background prediction. 



= 7.99 ± 1.05 pb 

Observed and expected ^-tagged jets after selection for 
different jet multiplicities are shown in Fig.jSj along with 
the expected contribution of inclusive tt signal, normal- 
ized to the measured cross section. 



IX. CONCLUSIONS 

We presented a measurement of the tt production cross 
section in a final state with large missing transverse en- 
ergy and multiple jets. We explicitly vetoed well iden- 
tified high-px electrons or muons from W boson decay, 
and rejected events with low missing transverse energy. 



to avoid overlaps with other cross section measurements 
performed by the CDF Collaboration. Using an opti- 
mized neural network based kinematic selection and a 
b jet identification technique on a sample of data cor- 
responding to an integrated luminosity of 2.2 fb^^, we 
obtain a production cross section value of 7.99 ± 1.05 pb, 
in good agreement with the reference theoretical value of 



7.46 



+0.66 
-0.80 



pb |lj and with other next to next to leading or- 



der calculations [28| for a top mass of 172.5 GeV/c^. The 
agreement with the most recent experimental determina- 
tions is also good. Given the high precision of the result, 
this independent cross section measurement will have a 
significant impact on a future CDF combined value. The 
same + j'^ts selection can be used for top quark mass 
measurement to obtain a result statistically uncorrelated 
with those of other methods. 
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